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Abstract 

We develop a one-parameter family of static baby Skyrme models that do not require a potential 
term to admit topological solitons. This is a novel property as the standard baby Skyrme model 
must contain a potential term in order to have stable soliton solutions, though the Skyrme model 
does not require this. Our new models satisfy an energy bound that is linear in terms of the 
topological charge and can be saturated in an extreme limit. They also satisfy a virial theorem 
that is shared by the Skyrme model. We calculate the solitons of our new models numerically 
and observe that their form depends significantly on the choice of parameter. In one extreme, 
we find compactons while at the other there is a scale invariant model in which solitons can be 
obtained exactly as solutions to a Bogomolny equation. We provide an initial investigation into 
these solitons and compare them with the baby Skyrmions of other models. 
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I. INTRODUCTION 


The baby Skyrme model mm is a nonlinear field theory admitting topological solitons 
known as baby Skyrmions. It is often studied as a (2+l)-dimensional analogue of the Skyrme 
model |3] for nuclear physics though is itself an interesting physical model with applications 
in condensed matter physics HHS]. In the (3+l)-dimensional Skyrme model, the topological 
solitons are called Skyrmions and can be used to model atomic nuclei with their topological 
charge, an integer B, giving the baryon number. As a lower-dimensional version of this 
model, the baby Skyrme model has been used to investigate a variety of difficult problems 
in the Skyrme theory including Skyrmion scattering HIM and the effect of isorotation on 
Skyrmion solutions [ME]. 

A key difference in the models arises when we consider their necessary components. The 
baby Skyrme model is an 0(3)-sigma model extended by the addition of a term quartic 
in derivatives called the Skyrme term and a symmetry breaking potential term. The com¬ 
bination of Skyrme term and potential gives a scale to the model and enables it to evade 
Derrick’s theorem [TB] for scalar field theories in two space dimensions. By contrast, in the 
full Skyrme model the combination of the Skyrme and sigma terms is sufficient to evade 
Derrick’s theorem. This provides one motivation for our paper — we wish to design a static 
baby Skyrme model that does not require a potential term to have topological solitons. One 
approach is to apply a noncommutative deformation to the baby Skyrme model instead of 
including a potential term We apply a very different method. 

Before discussing our approach in detail, we briefly review different variants of the baby 
Skyrme model. One way to create new models has been through the use of different poten¬ 
tials [ME], and it has been found that the choice of potential has a dramatic effect on the 
solitons of the model. In particular, the appearance and structure of multisolitons depends 
strongly on the potential term. For some potentials, higher-charge solitons form chains [20] , 
for some rings [TO] and for others [OTl |22[ stable multisolitons may not exist at all. Models 
have also been designed in which the 0(3) symmetry is broken to the dihedral group Djv, 
and here multi-Skyrmions have been observed with crystalline or broken structures [231 - 125] . 

In addition to choosing a different potential term, it is possible to develop new baby 
Skyrme models by removing the sigma term. Models consisting of only the Skyrme term 
and a potential are sometimes called restricted or BPS baby Skyrme models [2SH2H]. De- 
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formations of BPS models [29] have also been investigated, for which a physical motivation 
is found in the (3+l)-dimensional Skyrme theory. One signihcant problem in applying the 
Skyrme model to nuclear physics is that the binding energies of Skyrmions are considerably 
larger than the experimental values. The BPS Skyrme model [3nH32| consists only of a sextic 
term and a potential term, and has been developed, along with its generalizations, as an 
attempt to obtain more realistic binding energies. Different approaches to obtain Skyrme 
models with low binding energies are by coupling vector mesons to the 0(3)-sigma model 
term and removing the Skyrme term [2Sl El] or by studying Skyrme solitons on curved 
backgrounds [23 ES] • 

Another attempt to address the problem of obtaining realistic classical binding energies in 
the Skyrme theory has been to create new Skyrme models through a novel choice of potential 
term [371 EH] • A family of models is obtained by a one-parameter family of potential functions 
interpolating between the standard Skyrme model and a model in which a topological energy 
bound can be saturated for \B\ = 1. An equivalent idea has been explored in the baby 
Skyrme model to obtain so-called “aloof’ baby Skyrmions [2H]- We have also been motivated 
by recent interest in topological energy bounds [37100]. When designing our models, we 
require that they satisfy a particular topological energy bound. We hnd that this has several 
useful consequences for our models. 

Our approach to designing new baby Skyrme models is entirely different from those 
outlined above. As we wish to design models which do not require a potential term to 
have topological solitons, we remove the potential and raise the sigma and Skyrme terms to 
some powers. Considering Derrick’s scaling argument and requiring that our models satisfy 
a topological energy bound results in a one-parameter family of baby Skyrme models. We 
hnd that the required powers of the Skyrme and sigma terms are fractional. This draws a 
natural comparison between our models and the Nicole [H] and AFZ [l2l 03] models, which 
were investigated numerically in Refs. [HJ 02] along with a set of conformally invariant 
Skyrme-Faddeev models obtained by taking linear combinations of the two. 

The outline of the paper is as follows. We open in Sec. [TTjwith a brief overview of the static 
baby Skyrme model, focusing on its well-known energy bound and the application of the 
Derrick scaling argument to this theory. In Sec. uni we present our new models, beginning 
with a general form for the static energy and then illustrating how the application of Derrick’s 
theorem and the requirement that our solitons satisfy a topological energy bound reduces 
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the number of parameters to one. 

In the remainder of the paper we investigate the solitons of our one-parameter family of 


models. In Sec. IV we present our numerical results, hrst considering axially symmetric 
solutions and then progressing to simulations of the full held theory. We discuss the effect of 
the model parameter on the solitons and compare them to those found in existing models. 
We end by summarizing our results and rehecting upon open questions and opportunities 
for further investigation. 


II. THE BABY SKYRME MODEL 

The static energy functional of the baby Skyrme model is given by 

Ebs = j^ \ ^ + ^( 0 )^ ( 1 ) 

where the held 0 : —)■ 5^ is a three-component vector 0 = { 4 >i, 4 ’ 2 , 4 > 3 ) of unit length. 

The hrst term in Q is the 0(3)-sigma model term and is extended by the addition of a 
term quartic in derivatives, called the Skyrme term, and a potential term V (0) to allow the 
existence of stable topological soliton solutions. 

To ensure that solutions have hnite energy, the boundary condition 

lim 0= (0,0,1), (2) 

|x|^oo 

is imposed, assuming that (0,0,1) is a minimum of the potential V. This enables a one-point 
compactihcation U {cxd} = S^, and thus we can consider 0 as a map 0 : —)■ 5^. We 

can label the maps 0 by an integer B G vr 2 (S'^) = Z, called the topological charge. This is 
the winding number of the map, given by 

B = f 0 • (^10 X 020) d^x, (3) 

and is sometimes called the baryon number for comparison with the Skyrme model. The 
topological solitons of this theory are held conhgurations which minimize the energy ([^ in 
a given topological sector B. They are called baby Skyrmions. 

A lower bound on the energy of a solution with charge B in the baby Skyrme model is 
given by 

EBs>^7r\B\. (4) 
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This is a bound on the sigma term alone, obtained by a completing the square argument, 
and is never saturated by baby Skyrmions. 

When deriving energy bounds, it will be convenient for us to rewrite the static energy 
Q using its geometrical interpretation |16]. Dehne the symmetric, positive definite 2x2 
matrix D by 




(5) 


and let denote the eigenvalues of the strain tensor iA, where i = 1, 2. Then we can express 
the baby Skyrme energy functional Q in terms of the non-negative eigenvalues of D as 

Ebs = + A^) + 1a;A^ + V(4,)y^x, (6) 

and the topological charge can be expressed as 

B = -^ f X 1 X 2 d^x. (7) 

Using the energy (|^, we can obtain the well-known lower energy bound Q on the sigma 
term by completing the square as 


1 

2 


(A^ -I- A2) d^x 


1 

2 


(Ai =F X2f d^x ± 2 



> 


A1A2 d'^x 


4n\B\. 


( 8 ) 


The approach given above for deriving topological energy bounds is similar to those given 


when we derive energy bounds for our new baby Skyrme models. While Q is a well-known 
topological energy bound for the baby Skyrme model, recently tighter bounds have been 
obtained by also taking into account energy contributions from the Skyrme term and the 
potential term [271110]. 

The inclusion of a potential term in the baby Skyrme model is important as it allows 
the model to evade Derrick’s theorem [H] and thus have topological soliton solutions. This 
theorem rules out the existence of topological solitons in flat space scalar field theories by the 
requirement that a stationary point of the energy must also be stationary against rescaling. 
Therefore, if the energy of the theory after applying the spatial rescaling x h-/ rx, which we 
denote by e(/i), has no stationary point, then there can be no static finite energy solutions 
except the vacuum. 


in recent papers on the subject [371 HQ]. We will apply this method again in Sec. [Ill B 


5 





We apply this argument to the baby Skyrme model. Under the rescaling x h->■ px, the 
static energy Q becomes 

= E 2 + E 4 + fl '^Eq, (9) 

where we use E 2 , E 4 and Eq to denote the sigma term, Skyrme term and potential term, 
respectively. As a result of Derrick’s theorem, we observe that the combination of a potential 
term and the Skyrme term allows the existence of topological solitons. 

We can also derive a virial theorem satished by the baby Skyrme model by taking 
and setting this to zero, to hnd 

E 4 = Eq. (10) 

Contrast this with the results of applying the scaling argument to the Skyrme model, 
which has static energy 

Es = d^x, ( 11 ) 

where m is related to the pion mass, the pion helds are written as U : —)■ SU(2) and 

Ri = idiU)W. 

In this case, applying the rescaling produces 

es(h) = —E 2 + fiE^ H— ^Eq, (12) 

so the potential term is unnecessary to evade Derrick’s theorem. If we consider the static 
energy of the Skyrme model with no potential term, we can further obtain the virial theorem 


E2 — E4. 


(13) 


We have seen that the potential term is a necessary component of the baby Skyrme 
model if there are to exist topological soliton solutions. However, the same is not true of the 
Skyrme model, in which Skyrmions can exist without the presence of a potential term. This 
difference between the two theories motivates us to investigate the design of baby Skyrme 
models that do not include a potential term but still have soliton solutions. 

A variety of different functions V (0) have been investigated as the potential term in the 
baby Skyrme model [TBHTU] . Particular examples are 


U(0) 


Vi = m^(l — ^ 3 ) (“old” potential) 

V 2 = Tn?{l — 0 |) (“new” potential) 

V 3 = m^(l — ^ 3 )^ (holomorphic potential). 


(14) 


\ 
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The choice of potential has a strong effect on the appearance and structure of multisoliton 
solutions. For example, in the old baby Skyrme model Vi, higher-charge baby Skyrmions 
form chains [20] , in the new baby Skyrme model V 2 rings are minima, and in the holomorphic 
model V 3 no multisolitons exist. Recent work |39] has explored combining the old potential Vi 
with the holomorphic potential V 3 to obtain weakly bound multisolitons. We are interested 
in what the structure of baby Skyrmions would be without a potential to govern them. 


III. BABY SKYRME MODELS WITHOUT A POTENTIAL 

We propose a range of new baby Skyrme models that do not require a potential term 
to evade Derrick’s theorem. To achieve this, we raise the sigma and Skyrme terms to the 
power a and /9, respectively, and determine the range of acceptable values for these powers 
to ensure stability with respect to rescaling. As a starting point for the new static energy, 
take 

E = [ (ci{di(j) ■ dicj))"" + C2{\di(j) X (15) 

Jr2 V / 

where ci, C 2 are positive real coupling constants, and a, (3 are real constants. 


A. Derrick’s scaling argument 


To determine suitable values for a and /3, we apply the rescaling x ^ nx to the static 


energy (15) and consider the results of Derrick’s theorem. This leads to the energy 


e(/i) = 


(16) 


There are three cases in which our model can evade Derrick’s theorem: 

(i) a < 1 and f3 > 0.5, 

(ii) a > 1 and (3 < 0.5, (17) 

(hi) a = 1 and {3 = 0.5, 

with case (iii) providing a scale invariant model. We only consider cases (i) and (iii) because 

solutions in the models of case (ii) would either be compact or not have hnite energy; see 
Appendix for a detailed discussion. 
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There is also a virial theorem satished by our models. Take ^\^=i and set this equal to 
zero. For case (hi) this is automatically satisfied; otherwise, we have 

(2a - 2)^2 + (4/3 - 2)^4 = 0. (18) 

The resulting virial theorem is 

E 2 = E,. (19) 

a — 1 

Recall that the Skyrme model without a potential term satisfies the virial theorem 
E 2 = E 4 . Our models also satisfy this virial theorem when 

/S = 1 - f • (20) 

This selection of models includes one in which the static energy ( [l5| produces the same 
function 65 (/i) under rescaling as that for the Skyrme model ( [IT| ) without a potential term. 
In this case the parameters are a = 0.5 and /3 = 0.75. 


B. Energy bounds 


We have seen in Sec. that the baby Skyrme model ([^ satisfies a linear bound (|^ 
on the energy of its solutions in terms of the number of solitons. This is a useful property 
shared by many soliton models m- Therefore, we require our new baby Skyrme models to 
satisfy such a lower bound on the energy. In the following, we use this condition to fix the 


parameter /3 in (15) and further restrict the family of models that we consider. 


Defining the matrix D by Dij = dicp ■ djcj) as before, rewrite the energy (15) as 


E = 


Ci(A4 + A 2 )" + C2(2A4A2)^ ) d^X, 


( 21 ) 


where A^, denote the eigenvalues of D. 

To obtain a lower bound on the energy, we use the following special case of the inequality 
of the arithmetic and geometric means: for a, b non-negative, 

Qj b 


> y/ab, 


( 22 ) 


with equality if and only ii a = b. 






We obtain a lower bound on the energy by first applying inequality (|22|) twice to find 

E = 2 I 


-ci(Ai + A2)“ + ■^C2i2X\X\Y^(^x, 


>2 ^^(Xl + Xl)^^(2XlXl)^d^x, 

7r2 

> 2^+2i/ciC2 f 2t|AiA2|^|AiA2|^d^a;, 


(23) 


a+0 


a + 20 


= 2 ^+ 2 ^CiC2 / IA1A2I 2 d^x. 

JR2 

Then to ensure that this energy bound is linear in terms of the topological charge B, the 
required value of 0 is 


0 = 1 


a 


The resulting topological energy bound is 

E > 22 + 4yciC2 A'k\B\. 


(24) 


(25) 


Note that this choice of 0 was also found in Sec. Ill A by requiring that the virial theorem 


(19) be simply E 2 = Ei\ the virial theorem of the Skyrme model without a potential 


term 0. 

In case (hi), there is an alternative energy bound derived by a standard completing the 
square argument. The static energy for this model is given in terms of the eigenvalues A^, A^ 


as 


-^ = / {ci{X\ +Xl) + C2J~2>^^d^x. 

JR2 V V / 

By completing the square, we hnd 

E = [ ^cidAil — |A2|)^ + 2 C 1 IA 1 A 2 I + \/2c2|AiA2|')d^a;, 

JR2 ^ 2 

= [ ci(|Ai| - |A 2 |)^d^a; + (2ci + a/2c2) [ IA 1 A 2 I d^x, 

JR2 Jr2 

> 47r(2ci + V2c2)\B\. 

So an alternative bound in case (iii) is given by 

E > 47r(2ci + V2c2)\B\. 

This bound is saturated for solutions of the Bogomolny equation. 


(26) 


(27) 


(28) 


|Ai| — IA2I, 


(29) 
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which leads to the following system of equations: 

di(j) ■ di(j) = d 2 (j) ■ 024', (30a) 

di4 ■ O 24 = 0. (30b) 

Let z denote the complex coordinate z = x + iy in the spatial plane and R the Riemann 


sphere coordinate on the target S^. We can write solutions of Eqs. (30) in each topological 
sector B in terms of rational maps R{z) as 

1 


RR, —i{R — R), |-R| — 1 


(31) 


1 + 

Here R{z) = p{z)/q{z) is a ratio of two polynomials p{z) and q{z) with no common factors, 
and B = max{deg(p), deg(g)}. To satisfy the boundary condition lim|a;|_^oo 4 = (0,0,1) at 
infinity, we require R{oo) = 00 . One important case is the axially symmetric rational map 


R{z) = z^. 


(32) 


Thus, we can find exact solutions for baby Skyrmions of any charge B in this model. 

It still remains to set the values of the coupling constants ci and C 2 - In this paper, we 
choose 


_ 3 + Q: ^ 

Cl = 2 2 and C 2 = 22 ci, 

and so obtain the final form of the static energy for our models as 

E = 2~^ [ ({di4-0^4^+ 2{\di4x d24\‘^y~^)d‘^x. (34) 

JR2 ^ ^ 

This choice of coupling constants has three useful consequences. First, the choice of C 2 


(33) 


ensures that the bounds (25) and (28) in case (iii) coincide. It also causes the virial theorem 
E 2 = E 4 to be satisfied in case (iii), as we now have ci = \ and C 2 = and thus, 

E 2 = 2ci ■ 47r|i?| = 27r|R| = ■\/2c2 ■ 47r|i?| = F 4 . (35) 

Due to the scale invariance of this model, it is not necessary that the virial theorem be 


satisfied in this case. It is only due to the choice of constants (33) that the virial theorem 


(13) holds here. 

Finally, for any choice of a, this combination of ci and C 2 sets the topological energy 


bound (25) to be 


E > 4:7i\B\, 


(36) 
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with the bound saturated when a = 1. This is the well-known energy bound on the sigma 
term of the baby Skyrme model. 


IV. NUMERICAL RESULTS 


In this section, we calculate axially symmetric baby Skyrmion solutions for parameter 
a G [0.5,1] and with topological charges i? = 1 — 3, 10. We minimize the energy functional 


(34) for rotationally symmetric Skyrme conhgurations using two very different numerical 


approaches: ID gradient flow and Newton’s method for nonlinear systems. Finally, we 
perfom 2D energy minimization simulations for a selection of our models and verify that the 
minimal energy solutions agree with those found when imposing axial symmetry. Skyrmion 
chain solutions are found to be of higher energy. 

Note that the baby Skyrmion solutions for a = 0.5 are discussed in a separate subsection 
as the solitons become compactons. These are solitons with compact support, taking vacuum 
values everywhere outside some hnite region of space. Compact solitons have been studied 
before in the Skyrme-Faddeev model in the inhnite mass limit |18], and in massive baby 
Skyrme models [271 uni EQ]. Compactons are numerically challenging and require a careful 
adjustment of our numerical methods. 


A. Axial baby Skyrme solutions 


To hnd axially symmetric soliton solutions of the equations of motion, we use the ansatz 


0 = (sin/cos(56*), sin/sin(56*), cos/). 


(37) 


where r, 9 are the usual polar coordinates, / = /(r) is a radial prohle function, and B is 


the topological charge of the configuration. Substituting (|37|) into (|34|) yields the energy 

,sinV 


_ 3+a 

E = 2“^ ■ 27r 


/'^ + 5 " 


+ 2 f^B 


/2D2sin / 


2 rx i-f 


rdr. 


(38) 


which depends only on the radial coordinate r. Here prime denotes differentiation with 
respect to the radial coordinate r. By the principle of symmetric criticality, solutions of the 


Euler-Lagrange equation for the simplihed energy (38) will also solve the equations of motion 


for the original energy (34). In the following, we minimize the energy (38) by solving the 
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FIG. 1; Energy as a function of model parameter a for axially symmetric configurations with 
topological charges i? = 1 — 3, 10. (a) Total energy E. (b) Binding energy per soliton /\E/AttB. 


Euler-Lagrange equation (Al) in Appendix [^subject to the boundary conditions /(O) = vr 
and f\oo) = 0 in two ways: through the use of a ID flow method and also using Newton’s 
method for nonlinear systems [SI] with grid spacing Ar = 10“^ over the interval 0 < r < 20. 


For both methods, we monitor the topological charge and check the virial theorem (19) at 
each iteration step. 

Figure |l(a)| displays the total energy F as a function of model parameter a for axially 
symmetric baby Skyrmions with topological charges i? = 1 — 3, 10. The topological energy 
bound is indicated by a dashed line. All energy values are given in units of AnB. The energy 
values are furthest from the bound towards a = 0.5, but draw closer as a increases until the 
bound is hnally saturated in the a = 1 baby Skyrme model. As the charge increases, the 
difference between the energy value and the bound grows smaller. So the bound tightens for 
higher-charge solutions of a given model. Note that the energy difference between subsequent 
charges also decreases drastically. As a limiting case, we include in Fig. [Tjthe energy values 
for charge 10. In Fig. |l(b) we display the binding energy as a function of a, where the 
binding energy per soliton is given by 


AE _ Eb 
B ^ B ' 


(39) 


with El denoting the energy of the charge one solution and Eb denoting the energy of the 
charge B solution. The binding energy per soliton is the energy required to split a charge B 
baby Skyrmion into B charge one Skyrmions divided by the total number of solitons. The 
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a=0.51 



a=0.6 
a=0.7 
a=0.8 
a=0.9 
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FIG. 2'. B = 1 profile functions /(r) for model parameter 0.5 < a < 1. 

binding energy is fonnd to increase with the topological charge. 

Another featnre of the solntions which changes dramatically as a increases is illnstrated 
in Fig. Here we compare charge one profile fnnctions /(r) in a selection of the models 
ranging from a = 0.51 to a = 0.9 with the exact solntion of the Bogomolny eqnation for 
a = 1. The nnmerical profile fnnctions were calculated using the Newton method over the 
interval 0 < r < 20. For a = 0.51, the profile function is tightly concentrated between r = 0 
and r = 2.2. As a increases, the profile functions start to spread out. By a = 1, the profile 
function is less localized and approaches the vacuum gradually. 

A more detailed examination of the approach to the vacuum of the profile functions is 
given in Appendix In this appendix, we linearize the equation of motion as r —)■ cx) 
and obtain solutions that describe the profile functions as they approach zero. The profile 
functions exhibit a power law behavior, /(r) ~ for large r. As a tends to 0.5 the exponent 
A becomes increasingly negative, and the approach to the vacuum becomes steeper. At 
a = 0.5 the exponent diverges suggesting that solutions in this model are compactons. Near 
the origin, for all values of a the charge one prohle functions exhibit linear behavior. This is 
discussed in detail in Appendix where we linearize the equation of motion near the origin 
for any charge B. 
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B. Baby Skyrme solutions in the a = 0.5 model 


In this section, we present the results of 2D simulations for the compact charge one and 
two solitons obtained for model parameter a = 0.5. As a starting point for our 2D energy 
minimization routine we choose two different initial conditions: a rotationally symmetric 
conhguration created from the ID prohle function for a = 0.5 and the conhguration relaxed 
with a = 0.51. 


To hnd prohle functions in the a = 0.5 model we minimize the energy (38) over intervals 
r G [0, Rest] for various boundary points Rest surrounding the expected compacton radius. 
Solving the corresponding held equation over each interval is accomplished by Newton’s 
method for nonlinear systems with grid spacing Ar = 10“^ due to its increased speed over 
the gradient how method. We then seek the value of Rest that minimizes the energy. This 
enables us to decide upon a numerical energy value for the solution up to one decimal place 
of accuracy, with the virial theorem and topological charge also correct to one decimal place. 

By substituting the prohle functions obtained by this method into the axially symmetric 


ansatz (37), we build 2D conhgurations. These are implemented as initial conhgurations in 
a 2D relaxation method similar to that described in Ref. [10]. We evolve the equations of 


motion derived from (34) in a hctitious time t and include a damping term governed by the 
dissipation e. We periodically remove kinetic energy by setting 0 = 0 at all grid points. In 
the following, the hnite diherence approximations are second-order accurate in the spatial 
derivatives. The simulations are performed on a (401)^ grid with spacing Aa;=0.02 for the 
charge one soliton and Ax = 0.04 for the charge two soliton. In both cases the dissipation 
parameter e is set to 0.5. 


In Fig. 3(a) we plot the energy density of the resulting B = 1 conhguration. Figure 4(a) 


displays the energy density of a R = 2 conhguration obtained by the same method. The 
energy of both solutions is localized in a hnite region of space, and the steep approach to 
the vacuum is evident at the boundaries of the compactons. The energy values for these 
solutions agree with those of the corresponding prohle functions to one decimal place and 
are given in Table |T| 

The second method that we implement to hnd solutions in this model is to take a 2D 
conhguration with a = 0.51 as an initial condition in the 2D relaxation code. The same grid 


and spacing are chosen as for the previous initial conhguration. In Fig. 3(b) we compare 
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(a) 


(b) 


FIG. 3: Energy density for charge one baby Skyrmions with model parameter a = 0.5. (a) Surface 
plot of the energy density, (b) We compare slices through the energy density obtained when 
relaxing two different B = 1 initial conditions: a baby Skyrme configuration relaxed with a = 0.51 
(green line) and a rotationally symmetric configuration generated from an a = 0.5 profile function. 


charge one solutions obtained by both methods. We display slices along x = 0 through 


their energy density. The same comparison for charge two solitons is presented in Fig. 4(b) 


While the approach to the vacuum is not as steep for the second method, both methods 
generally agree well and describe the soliton’s energy to one decimal place. 


C. Higher charge solutions 


To verify our axially symmetric charge one and two solutions and to investigate solutions 
of higher charge, we implement a 2D numerical method. We apply the same relaxation 
method as in the previous section, but with a different grid. For models excluding a = 0.5, 
we use a (201)^ grid with spacing Ax = 0.2. To create initial conhgurations, we substitute 


our numerical prohle functions into the axial ansatz (37) to generate 2D conhgurations. We 


then take these conhgurations as initial conditions for our 2D energy minimization algorithm 
to hnd solitons of diherent a values. For example, an a = 0.8 axial solution is chosen as an 
initial conhguration to obtain the a = 0.7 solution. 

In Table [T] we give the energy values of our numerical simulations for a selection of a 
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FIG. 4: Energy density for charge two baby Skyrmions with model parameter a = 0.5. (a) Surface 
plot of the energy density, (b) We compare slices through the energy density obtained when 
relaxing two different B = 2 initial conditions: a baby Skyrme configuration relaxed with a = 0.51 
(green line) and a rotationally symmetric conhguration generated from an a = 0.5 profile function. 


values. All energy values are given in units of AttB, motivated by the energy bound of our 


models (36). We also present the binding energy for the 2D configurations, calculated using 


(39). For axial solutions the results of our 2D simulations agree to between two and three 


decimal places with the values obtained when minimizing (38). 

For higher charges, axial solutions remain the energetic minima, though other conhgu- 
rations have been obtained. In particular, we find chain configurations in our models by 
using three solitons in a line as an initial conhguration. The energy values for the chain 
conhgurations are also presented in Table and are denoted by a *. Their energy is higher 
than that of the axial conhgurations, and they do not satisfy the virial theorem. So these 
are local minima but not the global energy minimizers. 

In Fig. [^we compare the energies obtained by ID gradient how with those calculated by 
2D relaxation for baby Skyrmions with topological charges 5 = 1-3. We plot the energy 
for axially symmetric conhgurations with model parameter a = 0.51 - 1.0 and indicate all 
energy values computed by the full held simulations by points. As before, the topological 
energy bound is indicated by a black dashed line. The black points denote the energy of 
axially symmetric baby Skyrmions calculated by full held relaxation. They lie on top of the 
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FIG. 5: Total energy E for baby Skyrmions obtained with the axial ansatz (37) and for baby 


Skyrmions obtained with the 2D relaxation method. 


lines showing the energy for baby Skyrmions of the same topological charge obtained by 
ID gradient flow. The energy for charge three chain conhgurations is also included in the 
hgure as a series of red points connected by a dashed red line. For a = 0.8 the energy of 
the B = 3 chain is slightly higher than the B = 2 solution, but much lower than the B = 1 
solution. So, it is not energetically favorable for the chain to split up into a. B = 1 and 
B = 2 Skyrmion. 

The energy density for 5 = 1-3 baby Skyrmions with a = 0.6 - 0.9 is plotted in Fig. 
They are all axially symmetric, and the effect of increasing a on the solutions can be seen by 
comparing the graphs. This is most noticeable for the charge one solitons, where the energy 
density of the a = 0.6 solution is concentrated over a wide area with only a small tail. As 
a increases, the tail of the energy density becomes wider while the area in which the energy 
density is most concentrated decreases in width and increases in height. A similar effect 
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occurs for the charge two and three rings, which become thinner and taller as a increases. 

In Fig. we display the energy density of chain conhgurations for the same selection of 
models. The structure of the chain conhgurations changes signihcantly as a increases. For 
a = 0.6, the solitons are very close together but as a increases, the chain starts to pull apart. 
For a = 0.9, the chain almost splits into three separate solitons, though they remain close 
enough to deform each other. This may be explained by the approach to the Bogomolny 
solutions at a = 1. The attraction between solitons becomes weaker and weaker until they 
do not feel any attraction or repulsion at a = 1. Then the baby Skyrmions can be placed 
at arbitrary positions. In fact the solution space is the space of based rational maps. In the 
a = 1 model, the energy of three separate solitons is identical to that of a three-soliton ring 
solution. 

The chains observed at the lower a values most resemble those found in the baby Skyrme 
model [20], although the a = 0.6 chain in particular appears more squashed. The chains 
observed at higher a values are quite different. However, comparisons could be drawn 
between the a = 0.9 chain and the isospinning baby Skyrmions of ra in which chains are 
also seen to break up. 
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TABLE I: Energy values obtained from full field simulations /4ttB) when compared to ID 

gradient flow results /AttB). The binding energy per soliton AE/AttB is calculated by (39) 

using the numerical 2D energy results. The two B = 3 configurations given are the axial solution 
(37) and the chain configuration of Fig. which is denoted by 3*. 


a 

f3 

B 

/4ttB 

/47rB 

AE/47rB 

0.5 

0.75 

1 

1.2 

1.2 

0.0 

0.5 

0.75 

2 

1.1 

1.1 

0.1 

0.6 

0.7 

1 

1.188 

1.188 

0.0 

0.6 

0.7 

2 

1.092 

1.092 

0.096 

0.6 

0.7 

3 

1.068 

1.069 

0.119 

0.6 

0.7 

3* 

— 

1.081 

0.107 

0.7 

0.65 

1 

1.130 

1.130 

0.0 

0.7 

0.65 

2 

1.059 

1.058 

0.072 

0.7 

0.65 

3 

1.043 

1.042 

0.088 

0.7 

0.65 

3* 

— 

1.056 

0.074 

0.8 

0.6 

1 

1.068 

1.068 

0.0 

0.8 

0.6 

2 

1.028 

1.029 

0.039 

0.8 

0.6 

3 

1.020 

1.020 

0.048 

0.8 

0.6 

3* 

— 

1.036 

0.032 

0.9 

0.55 

1 

1.021 

1.020 

0.0 

0.9 

0.55 

2 

1.008 

1.007 

0.013 

0.9 

0.55 

3 

1.005 

1.005 

0.015 

0.9 

0.55 

3* 

— 

1.009 

0.011 
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FIG. 7: Contour plots of the energy density for B = 3 chain conhgurations. The value of a for 
each conhguration is indicated underneath its plot. 
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V. CONCLUSIONS 


We have developed a one-parameter family of baby Skyrme models that do not require 
a potential term to admit topological solitons. Starting with a general form for our models 
involving four parameters, we hxed three of these by specifying that our models should satisfy 
the topological energy bound E > 47r|i?|. Similarly to the Skyrme models described in [38] . 
we found that at one end of our parameter range there is a model in which this bound can 
be saturated. This model is scale invariant, and exact solutions to the Bogomolny equations 
can be obtained for any topological charge. Furthermore, our choice of parameters ensures 
that all of our models satisfy the same virial theorem, E 2 = E^, as the Skyrme model. In this 
way, we have designed a one-parameter family of baby Skyrme models without a potential 
whose scaling behavior better matches the Skyrme model and which even includes a baby 
Skyrme model scaling exactly like the Skyrme model at the extreme of our parameter range 
where a = 0.5. 

Our investigation into the solitons of our models showed that their form greatly depends 
upon the choice of the parameter a. Solitons in the a = 0.5 model are compactons. Both 
our numerical results for this model and linearization arguments support this conclusion. 
Previous examples of compactons in baby Skyrme models [271 SSI ED] depend on the choice of 
potential term used, and typically occur for particular parameter values in a one-parameter 
family of potential functions. There is no potential term in our models but the importance 
of parameter choice to the existence of compactons is similarly observed here. 

Solitons in the models with a > 0.5 were found to be less localized. We calculated 
solutions numerically for a selection of the models using three different methods. The energy 
of solutions to our models decreases as a increases, and higher-charge solutions are more 
tightly bound near the a = 0.5 end of the parameter range. As we approach the extreme 
of the parameter range in which the energy bound can be saturated, the binding energy of 
solutions decreases to zero. The energy minimizers are axially symmetric solutions, even for 
topological charge three. However, we also observed other higher-charge conhgurations with 
greater energy, in particular B = 3 chain conhgurations. Chain solutions for a G [0.6, 0.7] 
most closely resemble those observed in the old baby Skyrme model 1201 , while other chains 
have a very different appearance. In particular, as a increases our chain conhgurations begin 
to pull apart and become three almost separate solitons. 
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In this paper we have provided an initial study of our new baby Skyrme models but 
there is still further work that could be done. One interesting question is how the solutions 
behave for higher charges. Will the axially symmetric solutions always be the minimum 
energy configurations or will chain-like configurations play a more prominent role? An 
important challenge is to develop more accurate numerical methods to calculate compacton 
solutions, as this may have important applications beyond the area of topological solitons, 
for example in relation to fractional Laplacians |52] and fractional diffusion [53]. Although 
the unique aspect of our models is their lack of need for a potential term, the inclusion 
of a potential term in this setting offers many possibilities. Following [ID] and balancing 
the different terms with the potential would give rise to multi-parameter families of models 
which still obey a linear energy bound. The choice of potential has an important effect on 
the structure of solitons in the usual baby Skyrme model. Whether the same is true of our 
new models and generalizations thereof remains to be determined. 
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Appendix A: Linearizing as r —)> oo 

In this appendix, we linearize the equation of motion for an axially symmetric charge B 


soliton in the baby Skyrme model (15) for large r. As discussed in Secs. Ill A and IIIB, we 
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FIG. 8: Real solutions A+, A_ of Eq. (A2) as functions of a for 0 < a < 1 and B = 1. 


set /? = 1 — f, and the constants ci and C 2 are fixed in (33). 


The eqnation of motion for the prohle fnnction /(r) is given by 


4«(«-l)cir , 2/3(2/3-l)c2(2r^sinV)''(r+ ^sinV)'"“ 


r 


r + f sin^ / 


+ 


/ 


/2 


-|- ‘ 2 :OiC\ 


^ 2aci ^ 2a{a - l)ci ( - 2^ sin^ / + sin 2/) 


/'2 + f sinV 

^2 ■2r\Pfpf2 I B 2 2 


C 2 /? (2/3-1) {2p f sin^ fY{P + ^ sin^ /) f sin^ f + f sin 2/) 


/'sinV 


-aci^ sin2/ = 0. 


(Al) 


To linearize the symmetric eqnation of motion (Al) as r —>■ oo, we snbstitnte /(r) = r 


into (Al) and consider only the leading-order terms. For a < 1, this leads to the eqnation 
(2a - 1)A^ - 2(a - 1)X^ + 2B\a - 1)\‘^ - 2B\a - 1)A - B^ 


A 2 + 52 

which simplihes to the qnadratic 

(2a - 1)A2 - 2(a - 1)A - = 0. 

For a ^ ^ this can be solved for A, to hnd 

1 ± ^a 2 +2(52-1)0-52 + 1 


= 0 , 


A+ = 


a 


2a-I 


(A2) 


(A3) 


(A4) 
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For parameter value a = 1, most of the terms cancel and the solutions are 


A = ±B. 


(A5) 


At the other extreme, where a = 0.5, the quadratic term in (A3) cancels. The solution A_ 


becomes singular and A+ = . In Fig. [^we plot the two solutions A+, A_ of Eq. (A3) as 

a function of the model parameter a for 0 < a < 1 and topological charge B = 1. The root 
A+ is always positive and therefore we discard it. The root A_ has an asymptote at a = 0.5 


where Eq. (A2) is singular. For 0 < a < 0.5, the values of A_ are positive and therefore 


solutions for this range of a are either compact or do not have hnite energy. The interesting 
parameter range is 0.5 < a < 1 in which A_ < —1 and solutions have hnite energy. 


For a > 1, the leading-order terms in (Al) after substituting /(r) = are different, 
so that a different linearized equation is found in this parameter range. Here the equation 
becomes 


2(2 — a)(l — a)c2 , 
A 




2 \ 


(A 6 ) 


The only solutions are A = 1 and A = ±ii?. Thus for a > 1, solutions to the equation of 
motion would either have inhnite energy or be compact. 

Appendix B: Linearizing near r = 0 


We also linearize the equation near r = 0 to gain a greater understanding of the behavior 
of solutions in relation to their topological charge. Set 


/(r) = TT — ar'’'. 


(Bl) 


in Eq. (Al), where a is constant, and assume that 7 > 1. For small r'’', we can use the 


small angle approximation to replace trigonometric terms. By our assumption on 7 and the 
constraint 0.5 < a < 1, we hnd the leading-order equation 

,3/„, 1-7 , .A 1 ^D 2 7-1 , d 2 


2aaci ( 27 '^(a — 1) 


52 + 72 

which simplihes to the quadratic 


+ 7(1 - 7 ) - 27(0 - 1)5" 


52 _|_ ^2 


+ 52-7 =0, (B2) 


52 + 2{a - 1)7 - (2a - 1)72 = 0. 


(B3) 
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For a > 0.5, we can solve this to find the positive root 


—1 + a + + 2{B‘^ — l)a + 1 — 

2a-I 


(B4) 


Note that once again the a = 0.5 case mnst be considered separately. Here we obtain the 
leading-order eqnation 


B -7 = 0, 


(B6) 


SO we find 7 = when a = 0.5. This is also the limit of (B4) as a —)■ 0.5. 


At the other end of onr allowed range of a values, where a = 1, we observe that the 


expression (B4) simplifies to 


7 = \B\. 


(B6) 


Notice that for B = 1 and any choice of a > 0.5, we find 7 = 1 . So for all of onr models, 
the charge one profile fnnction has a linear behavior near the origin. For B > 1 we confirm 


that 7 (a) as given by (B4) does not have any tnrning points in the interval 0.5 < a < 1, 
with 7 ( 0 . 5 ) = B^ and 7 ( 1 ) = \B\. Hence 7 > 1 for any choice of H > 1. This justifies our 
earlier assumption on 7 . 
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